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Clifford  J.  Nolan  William  W.  Symes 

May  19,  1995 


Abstract 

We  consider  scattering  associated  to  the  reduced  scalar  wave  equation.  High  fre¬ 
quency  asymptotic  solutions  of  this  equation  leads  to  the  theory  of  geometrical  optics. 
In  this  theory  energy  is  transported  along  rays  (orthogonal  trajectories  to  wavefronts). 
However  this  theory  breaks  down  as  soon  as  the  ray  field  forms  an  envelope  called  a 
caustic.  This  signals  that  a  dramatic  change  in  nature  of  wave  propagation  occurs  in 
the  vicinity  of  a  caustic.  To  illustrate  this  change  of  character  we  study  an  experiment 
which  shows  that  reflected  waves  may  have  arbitrarily  high  energy  content  relative  to 
the  “size”  of  the  scatterer.  Moreover  a  theorem  is  proved  showing  that  this  unbounded 
behaviour  can  only  occur  when  a  caustic  develops. 


1  Introduction 

Our  motivation  for  this  study  is  the  scattering  of  acoustic  waves  in  the  subsurface  of  the 
earth.  However  the  results  are  completely  general,  applying  to  any  physical  phenomenon 
described  by  the  scalar  wave  equation.  We  shall  set  up  two  physical  models  to  demonstrate 
applicability  of  our  results.  First  we  describe  the  scattering  problem  in  geophysics  and 
then  an  analagous  problem  in  electromagnetism. 

For  simplicity  we  consider  the  earth  as  a  two  dimesnsional  halfspace  (although  the 
results  in  this  paper  can  be  generaized  to  any  dimension),  with  cartesian  coordinates 
(x,y)  occupied  by  a  fluid  which  has  constant  density  and  propagates  pressure  waves  in 
its  interior  at  a  local  speed  c(x,y).  The  Green’s  function  p(x,y,t)  associated  to  pressure 
waves  travelling  in  the  subsurface  of  the  earth  as  a  result  of  an  energy  source  with  time 
signature  /  G  L2(R)  being  fired  on  the  surface  at  location  x,  satisfies  the  acoustic  scalar 
wave  equation 

c~2(x,y)  ^y( x,y,t)~  A p(x,y,t)  =  f(t)6(x-xs).  (1.1) 

Similarly  in  the  electromagnetic  case  we  have  an  electromagnetic  field  propagating  in  a 
two  dimensional  half  space  as  above.  The  electromagnetic  field  obeys  Maxwell’s  equations. 
Let  D,  B,  E,  H,  J  respectively  denote  the  displacement  current,  the  magnetic  field,  the 
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electric  field,  the  magnetic  intensity,  and  the  current  density.  Then  Maxwell’s  equations 
are 


V  •  D  —  Pfree 


VB  =  0 


V 


x  E  =  - 


3B 

dt 


V 


X  B  =  J /ree 


In  addition  to  these  equations  we  have  the  constitutive  equations 

J  =  ctE 


(1.2) 


B  =  (poM) 


D  =  c0cE  (1-3) 

where,  e  is  the  relative  electric  permitivity  tensor  and  Co  is  its  value  in  free  space.  Similarly 
fi  is  the  magnetic  permeability  tensor  and  po  is  its  value  in  free  space. 

Now  we  make  some  simplifications  to  get  to  the  reduced  wave  equation  for  the  electro¬ 
magnetic  field. 

Assume  that  the  electromagnetic  field  is  studied  in  a  2  dimensional  half  space  which 
is  electrically  neutral  (Pfree  =  0),  non-magnetic  (p  =  po)>  and  non- conductive  (a  =  0). 
Coordinates  on  this  half  space  will  be  cartesian  and  a  typical  point  will  be  denoted  by 
(x,y)  and  in  future  we  will  often  write  x  in  boldface  for  ( x,y ). 


With  these  assumptions,  it  then  follows  [7]  from  Maxwell’s  equations  (1.2)  that  for  an 
electric  field  polarized  along  the  x-axis 


E(x,y,t)  =  (ip(x,y,t),0) 


(1.4) 


satisfies  the  reduced  wave  equation  equation 


c  2(x,y) 


d2rp 

dt 2  1 


x,y,t)~  At l>(x,y,t)  = 


where  0  <  c(x,  y) 


leo{x,y)no(x,y) 

*{x,y) 


(1.5) 


At  this  stage  we  start  our  general  study  of  the  scalar  wave  equation. 
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We  study  the  related  Green’s  function  for  the  reduced  wave  equation 

c-2(x,  v)  ^(x,«)  -  AG(x,l)  =  «(l)«(x  -  *.)  (1-6) 

corresponding  to  an  impulsive  source  at  location  (x»,0)  -  see  figure  (1.1)  for  the  corre¬ 
sponding  experiment.  This  enabales  us  to  find  if),  since  if)  is  just  the  right  hand  side  of 
(1.5),-  convolved  with  G. 


Figure  1.1:  The  experiment 


Assume  that  the  coefficient  c  has  been  split  into  two  pieces 

c(x)  =  c0(x)  +  6c(x) 


(1.7) 


where  c0  is  a  smooth  function  and  Sc  has  much  higher  frequency  content  than  cq. 


Denote  by  Go  the  solution  of  (1.6)  with  c  replaced  by  Co.  In  (1.6),  write  G  —  Go  +  SG 
to  define  6G.  Then  if  one  linearizes  equation  (1.6)  about  c0,G o  one  finds  that  SG  satisfies 
the  linearized  wave  equation 


1  d26G 
c2  dt2 


-  v  •  (ysG)  = 


2 Sc  d2G0 
eg  dt2  ' 


(1.8) 


Thus  Sc  appears  as  a  source  for  the  (scattered)  wavefield  SG.  Then  Go  +  SG  is  the 
Born  approximation  to  G. 

Our  ultimate  goal  is  to  solve  the  following  inverse  problem.  Given  that  we  record 
the  scattered  field  on  a  finite  portion  [0,  L]  of  the  line  y  =  0,  deduce  Sc.  Crucial  to  the 
solving  of  such  an  inverse  problem  is  a  good  understanding  of  the  nature  of  the  forward 
problem  i.e.  the  operator  which  maps  Sc  into  such  a  recording  on  y  =  0.  We  now  formally 
define  this  operator  and  proceed  to  investigate  its  character  in  the  later  sections.  We 
do  not  emphasize  the  inverse  problem  in  this  paper,  but  instead  make  reference  to  the 
implication  of  our  results  for  it. 

Definition 

( P6c)(x,t )  =  m(x,t)  • 6G(x,0,t )  (1.9) 
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where  m  €  Cq °(Rx  x  [0,T])  is  a  cutoff  function  which  is  used  to  reflect  the  fact  that  we 
only  record  the  scattered  field  for  a  finite  time  T  and  for  a  finite  portion  [0,  L ]  of  the  line 
y  =  0. 

Note  that  we  are  focussing  on  the  linear  behaviour  of  T  as  a  function  of  the  perturbation 
Sc  alone.  Thus,  tacitly  we  assume  that  we  have  estimated  cq  (  there  are  various  techniques 
available  for  this  )  in  order  to  get  an  expression  for  Go  (actually  an  asymptotic  expansion 
for  Go  will  be  used)  as  we  shall  later  see. 

The  Main  Result 


We  will  show  that  if  we  impose  some  mild  restrictions  on  the  geometry  of  rays  travelling 
in  the  medium  defined  by  the  coefficient  co,  then  we  have  that  for  a  certain  Hilbert  space 
H  that 

->  L2([0,I]x  [0,T])  (1.10) 

boundedly  iff  and  only  if  the  rays  in  Co  do  not  form  caustics.  This  result  will  be  stated 
more  precisely  as  a  theorem  in  section  4. 

Consequence 

Besides  the  obvious  physical  interpretation  of  the  result,  the  theorem  has  the  following 
consequence  for  the  inverse  problem.  One  needs  to  have  knowledge  of  the  spectral  prop¬ 
erties  of  T  for  an  efficient  inversion  algorithm.  The  theorem  shows  that  T  may  have  a 
large  spectral  radius  when  caustics  are  present  in  the  medium.  This  indicates  a  source  of 
trouble  for  inversion  algorithms  only  using  local  information,  for  example  Gradient  based 
algorithms  used  to  find  the  minimum  of  an  objective  function. 

Finally  for  reference  sake  we  give  a  list  of  the  Green’s  functions  in  various  dimensions 
in  a  homogeneous  medium  cq  =  constant. 


n  =  1:  G0(x,;x,t)  =  H(c0t  —  |x-  x,|) 
Co  H{cot  -  lx  -  x3|) 


n  =  2  :  G0(x5;x,t)  = 


2*  yjclt*  -  |X  -  X,| 


n 


3:  G0(xf;x,t) 


S(c0t  — Jxj— Xgj) 
47T  jx  -  x,| 


(1.11) 


2  The  progressing  wave  expansion 


As  promised  in  the  previous  section,  given  that  we  know  Co,  we  can  give  an  asymptotic 
expansion  of  the  Green’s  function  Go-  This  expansion  is  known  as  the  progressing  wave 
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expansion  (see  [3],  page  188).  The  description  naturally  leads  to  the  notion  of  a  caustic 
which  we  define  in  this  section  too. 

According  to  this  theory  we  have  an  expansion  summarized  by 

G0(x,;x,t)  =  a(x,;x)5(t-  r(x,;x))  +  Smoother  terms. 

Here  the  amplitude  a  satisfies  the  transport  equation 

V-(  a2(xJ;x)Vr(x4;x) )  =  0  (2.1) 

and  the  traveltime  r  satisfies  the  Eikonal  equation 

|Vr(x4;x)|2  =  (2.2) 

Also  we  included  xs  just  as  a  parameter  in  the  above  expressions. 

Now  we  will  discuss  the  meaning  of  the  above  symbols.  S(x,t)  is  the  most  singular 
term  appearing  in  the  Green’s  function  for  a  homogeneous  (c(x)  =  constant)  medium. 
For  example  in  one  dimension  S  is  given  by  (see  the  list  of  Green’s  functions  given  in  the 
last  section) 

S(x„x,t)  =  H(cot  -  |x-x,|) 

etc.  The  level  curves  of  the  traveltime  function  r(x*;x)  are  just  the  wavefronts  emanat¬ 
ing  from  the  source  location  x„.  Alternatively,  instead  of  dealing  with  wavefronts,  we 
can  consider  the  equivalent  concept  of  rays  which  are  the  orthogonal  trajectories  of  the 
expanding/contracting  wavefronts.  Then  r(x,;x)  =  r(x)  represents  the  amount  of  time 
that  a  ray  takes  to  travel  from  the  source  x,  to  x.  Similarly  along  the  rays  there  is  an 
associated  amplitude  a(x„;x)  of  the  signal.  This  represents  that  amplitude  of  the  signal 
at  x  which  departed  from  x,  along  a  ray  which  intersects  x. 

Of  course  implicit  in  all  of  the  above  description  of  r(x4,  x)  and  a(xa;  x)  is  the  fact  that 
we  assume  that  there  is  only  one  ray  leaving  x,  and  arriving  at  x;  i.e.  we  assume  these 
functions  are  single  valued.  This  is  usually  only  locally  true  for  even  mildly  inhomogeneous 
media.  When  this  fails  to  be  the  case  there  is  generally  a  locus  or  curve  to  which  there 
is  a  ray  tangent  at  every  point.  For  example  see  figure  (5.9)  for  a  distinctive  cusp  shaped 
curve  with  this  property.  Such  a  curve  is  called  a  caustic.  A  more  precise  mathematical 
description  of  a  caustic  is  given  by  reference  to  the  material  in  [4],  section  3.6.  For  the 
reader  who  wants  a  more  precise  description  of  a  caustic,  there  follows  a  summary. 

The  Eikonal  equation  is  first  order  nonlinear  partial  differential  equation  (p.d.e.)  and 
we  recall  how  local  solutions  to  such  equations  may  be  built  geometrically  via  the  method 
of  characteristics  as  follows. 

Consider  solving  the  following  first  order  nonlinear  (p.d.e.)  described  by 

/(x,f^)  =  °>  x  €  X; 
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tIxCx')  —  s(x ); 
Vt|a4(x')  =  h{x') 


(2.3) 


,  ,  ,  /  _  /  /  ri  \  are  appropriate  coordinates  on  the  hypersurface  M. 

for  r(x),  where  x  -  (in-iVi)  PP  p  Vvnprcnrface  M  If  for  all  points  x'  on 
rt  fipnote  the  normal  vector  field  on  the  hypersurlace  M-  ^  v  . 

Let  i/(x )  denote  tne  nor  then  M  ig  called  a  non  characteristic 

the  hypersurface,  we  have  that  /(x  ,  f^x  ))  /=  > 

hypersurface  for  the  p.d.e.  (2.3). 

The  method  ot  characteristics  say.  that  if  M  is  a  non  characteristic  for 

the  p  d.e  (2.3),  then  a  local  solntion  of  (2.3)  near  M  can  be  constructed  as  follows. 

Consider  a  neighbourhood  (in  *")of  Af  which  we  will  call  X.  We  will  define  a  solution 
of  (2.3)  in  an  open  subset  of  X. 

We  denote  by  T‘X  the  cotangent 

in  T-XwWch  we"  wffl 'afsume  is^sl"^™  trying  to  And  a  local  solution  to  equation 
2.3)  if  equivalent  .0  finding  the  sero  se,  of  the  so  called  Hamiltonian  nnc  m.  /(x,£) 

-  > - * 

[4]  for  a  justification  of  the  foUowing  prescription  for  a  solution. 

Associated  to  the  Hamiltonian  function  /(*,{)  is  a  vector  field  H,  on  T'X  given  in 
the  above  coordinates  by 


(2.4) 


Here  {  £lw>  }’  «™  *  *he  *“»-*  Spa“  *°  PX  * 

that  i^t)  •  /  =  0  verifying  that  /  is  constant  (equd  .0  zero)  along  the  Sow  o,  the 

vector  field  H j. 

bicharacteristic  curves  How  as  long  as  the  resulting  A  r™““  j™',™', lore  A  ha3  focal 
in  T’X,  then  we  will  have  the  following  situation.  /U  =  V  d  for  some 

coordinates  (%,{(*»,  »l>«e  «x)  is  a  smooth  function  given .by -  £ (  >  -  ^  , 

smooth  ^  ^  “ 

caustic. 
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What  is  interesting  to  us  is  what  happens  when  this  transversality  condition  mentioned 
in  the  last  paragraph  fails  to  hold.  This  says  that  the  canonical  projection  II  :  (x,  £)  — *•  x 
from  A  into  X  fails  to  be  a  diffeomorphism.  i.e.  when  projecting  from  a  neighbourhood  of 
such  a  location,  then  many  bicharacteristic  curves  get  “concentrated”  in  X  upon  projec¬ 
tion.  The  projection  of  a  bicharacteristic  curve  into  X  is  known  as  a  characteristic  curve 
or  simply  a  ray.  The  rays  mentioned  before  are  exactly  this  concept  applied  to  the  Eikonal 
equation.  If  there  is  a  subset  C  in  X  which  is  the  image  of  the  canonical  projection  from 
A  and  the  projection  is  not  a  diffeomorphism  here,  then  this  curve  X  is  an  envelope  of  rays 
by  definition.  Under  these  circumstances  X  is  known  as  a  caustic.  In  the  case  n  =  2  and 
X  being  a  curve  (not  necessarily  smooth)  this  curve  is  known  as  a  caustic  curve. 

We  give  the  following  schematic  (since  we  cannot  give  a  4  dimensional  picture)  figure 
as  an  example  of  the  above  ideas. 


3  A  motivational  experiment 


We  will  now  give  an  argument  for  the  main  result  of  the  paper.  The  motivation  consists 
of  two  parts.  First  we  give  a  methematical  reason  which  strongly  sugests  that  the  result 
may  be  true,  and  then  we  give  a  heuristic/physical  motivation. 


3.1  The  mathematical  plausability 

We  shall  soon  see  that  the  forward  operator  T  is  a  Fourier  Integral  Operator  [4]  and  we 
shall  quote  a  result  about  such  operators  which  is  the  basis  of  the  argument  here.  We 
state  the  result  here  for  convenience. 

Theorem  3.1  (See  -  [4],  page  154)  Let  C  C  T*X  X  T*Y  -  {0}  be  a  homogeneous 
canonical  relation  such  that 
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1. 

{(0  *  (  0  )  }  ^  ^(*,(,1/,*)  (^0> 

2. 

{(S)  -  (?)}  *  (C) 

3.  The  differential  of  the  projection  C  — +  T*X ,  (reap.  C  — *  T*Y )  has  rank  >  dim  X  +  k 
(resp.  >  dim  Y  4-  k).  Here  k  £  Z. 


Then  a  Fourier  Integral  Operator  A  £  I'T'(X,Y\C)  is  continous: 


H 


Comp 


(Y) 


„3-m-±(dim  X+dim  Y-2k)  ^  „ 

nLoc  ,  Vs,  m  £  H. 


(3.1) 


In  our  case  it  will  turn  out  that  the  forward  operator  T  £  /^([O,  L\  x  [0,T],X,  ;A) 
for  X  equal  to  a  half  space  representing  the  earth’s  subsurface  and  A  is  generated  by  a 
family  of  bicharacteristics  associated  to  the  wave  equation  -  see  [12].  It  so  happens  that 
when  caustics  develop  in  the  ray  field,  that  the  number  k  mentioned  above  drops  from 
k  =  dim  X  =  dim  Y  to  some  smaller  integer.  When  this  happens,  the  index  on  the 
righhand  side  Sobolev  space  in  (3.1)  is  strictly  less  than  s  —  m.  In  our  case  this  means 
that  if  a  caustic  is  present  then  the  theorem  is  telling  us  that  we  cannot  be  sure  that  T  is 
continuous  from  H*(X)  to  L2([0,  L\  x  [0,T]).  We  can  therefore  reasonably  hope  to  prove 
that  T  is  not  continuous  in  the  presence  of  caustics  and  consequently  that  the  estimate 
in  the  theorem  is  sharp. 


3.2  The  physical  plausability 

This  argument  is  based  on  the  following  experiment  which  was  first  suggested  in  [2].  It 
was  also  this  work  which  first  suggested  that  a  theorem  like  the  one  mentioned  in  the  last 
section  may  be  true.  In  the  former  thesis  [2],  formal  arguments  are  presented  in  support 
of  such  a  theorem.  The  ideas  presented  in  this  section  are  also  needed  for  the  next  section. 

Here  we  consider  performing  a  sequence  of  separate  experiments  labelled  by  an  index 
j,  whereby  a  reflectivity  ( 6cj )  is  introduced  to  reflect  the  incident  bundle  of  rays  back  to 
the  surface  of  the  earth  in  the  form  of  a  plane  wave  (see  figure  (3.1)).  Before  proceeding 
any  further  we  define  more  precisely  how  the  “reflectivity”  sequence  Scj  is  to  be  chosen. 
We  will  restrict  our  attention  to  the  reflectivity  sequence  of  the  form 

*Cj(x)  =  Xj(x)  eik>  x  (3.2) 
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where  Xj  18  a  test  function. 

Next  we  have  to  describe  the  functions  Xj  while  we  reserve  the  choice  of  sequence  of 
vectors  kj  until  later.  In  order  to  describe  the  supports  of  the  reflectivities  Xj,  let  us  denote 
by  A  the  curve  prescribing  a  caustic  that  develops  as  a  result  of  tracing  rays  through  the 
medium  from  our  point  source  x,  -  see  figure  (3.1).  Also  let  us  denote  by  T  the  interior 
of  the  closed  ray  tube  generated  by  the  rays  leaving  the  the  source  point  x,  and  traveling 
until  the  caustic  curve  A  is  touched  by  them.  We  will  refer  to  T  as  the  “precaustic  region  . 
In  the  next  section,  we  will  give  a  more  mathematical  description  of  the  precaustic  region. 
For  the  general  element  in  our  sequence  we  set  Xj  e9ual  to  t^e  bump  function  supported 
on  a  disk  of  radius  Vj  and  with  centre  (xj,yj)  i.e. 

f  if  y/(x  -  Xj)2  +  (y-  yj)2  <  vj 


A  |  0  ,  if  J(x  -  Xj)2  +  (y  -  yj)2  > 

wViptp  - 

4  =  -  V(*-*i)a +  (»-«)’■ 

The  quantities  Vj  are  chosen  to  be  a  strictly  monotonic  decreasing  positive  sequence,  while 
Xj ,  yj  are  chosen  such  that  the  distance  between  the  centre  point  (xj,y j)  and  the  caustic  A 
is  also  strictly  decreasing.  We  will  require  that  these  disks  do  not  touch  the  point  x,  and 
are  also  contained  in  T.  Furthermore  we  will  refer  to  the  interior  of  Supp  Xj  as  fy-  Thls 
situation  is  depicted  in  the  figure  (3.1).  Note  that  the  same  bundle  of  rays  is  intersected 
(“trapped”)  by  the  reflectivity  8cj  in  each  experiment  (i.e.  for  all  j).  Moreover  the  radius 
of  the  disk  flj  — ♦  0  as  its  centre  approaches  the  caustic  (i.e.  for  all  j). 

The  fact  that  it  is  always  possible  to  “trap”  the  same  bundle  of  rays  and  reflect  them 
into  a  plane  wave  says  something  about  the  shape  of  the  reflector  in  order  that  Snell’s  law 
be  satisfied.  That  this  experimental  set  up  is  possible  (locally)  is  simple  a  consequence  of 
the  Implicit  Function  Theorem  (see  the  appendix  for  a  proof). 

The  experiment  also  suggests  that  the  reflectivity  sequence  reflects  arbitrarily  large 
amounts  of  energy  relative  to  the  “size”  of  the  reflector  support.  Recall  that  the  wave 
energy  in  a  region  is  proportional  to  the  L2  norm  over  said  region  of  the  amplitude  of  the 
signal.  In  the  precaustic  region  close  to  a  caustic  the  raytube  T  collapses  (its  diameter 
tends  to  zero)  and  Geometrical  Optics  (assumed  valid  here)  predicts  that  the  amplitudes 
“blow  up”.  The  reader  can  verify  this  immediately  by  integrating  equation  (2.1)  over 
the  raytube,  applying  the  divergence  theorem  and  noting  that  Vr  is  parallel  to  the  rays 
generating  the  raytube.  This  predicts  an  arbitrarily  large  energy  density  as  we  approach 
the  caustic  and  the  raytube  collapses.  However  there  is  a  limitation  to  the  above  argument 
as  one  might  have  suspected.  The  above  argument  relied  on  geometric  optics  (ray  theory) 
which  is  only  valid,  provided  we  are  on  the  order  of  one  wave  length  of  the  incident  field 
away  from  the  caustic.  There  exists  a  more  complicated  uniform  expansion  (see  [16], [5]) 
in  the  neighbourhood  of  a  caustic. 
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Our  source  function  S(t )  in  equation  (1.6)  contains  all  frequencies,  and  hence  arbi¬ 
trarily  small  wavelengths.  Thus  we  may  expect  geometrical  optics  to  be  valid  for  these 
components  arbitrarily  close  to  the  caustic  A. 

Based  on  this  experiment  one  may  then  suspect  that  the  forward  map  T  is  an  un¬ 
bounded  map,  i.e.  that  we  can  find  a  reflectivity  sequence  such  that  the  ratio 

WFtcjW. 

\\Scib 

as  the  experiment  index  j  — *  oo.  The  subscripts  a  and  (3  indicate  two  different  norms;  the 
choice  of  which  will  become  clear  in  the  next  section.  Our  task  is  to  show  by  example, 
unboundedness  of  the  latter  ratio  and  hence  unboundedness  of  the  forward  map.  This  is 
the  main  result  of  the  paper. 


Figure  3.1:  A  Sequence  of  reflectivities 


4  Unboundedness  of  the  forward  map 

In  this  section  we  give  a  rigoroous  proof  of  the  main  result  of  the  paper.  Let  us  first 
establish  the  hypotheses  and  notation  for  this  section.  Let  Cl  be  some  open  subset  of  Rn . 
Then  define  JT*(fi)  as  the  closure  test  functions  Co°(ft)  in  H *. 

Throughout  this  section,  a  standing  hypothesis  will  be  made.  It  will  be  referred  to  as 
“the  simple  ray  geometry  hypothesis”  (SRG),  and  we  will  now  describe  it. 

Recall  from  section  2,  that  associated  to  solving  the  Eikonal  equation  (2.2)  in  Rn , 
there  is  a  2n  system  of  ordinary  differential  equations  (O.D.E.’s).  This  system  is  known  as 
Hamilton’s  Equations  and  solving  it  amounts  to  computing  a  Hamiltonian  flow  associated 
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to  the  function  on  /(x,  £)  =  j£|2  -  c-2(x).  The  system  of  autonomous  O.D.E.’s  follows 

x=  V*/(x(t),f(t)) 

£=  -vx/(x(t),f(t)) 

where  •  =  -7-.  (4.1) 

at 

Let  9  €  S1  and  denote  by  (  Xs(t),  £s(t)  )  the  solution  of  equations  (4.1),  which  have  initial 
conditions  (  x4,  ).  Note  that  t  is  just  the  ordinary  differential  equation  parameter. 

Define  the  following  sets 

0(x4)  =  |  x  :  x  =  Xg(t ),  for  some  9  6  Sn~l,  and  t  €  [0,T]  } , 

fl°(x4,  t)  =  {  x  €  0(x4)  :  if  x  =  Xg(t),  then  for  0  <  t  <T,  Xg{t) 
lies  only  on  the  ray  Xg  and  r(x4,x)  <  t  } . 

The  boundary  points  of  f belong  to  caustics,  and  Sl°(x3,T)  defines  more  precisely 
what  we  have  been  calling  the  precaustic  region. 

Now  we  state  the  SRG  hypothesis: 

1.  All  reflectivity  fields  are  contained  in  fi°(x,,T). 

2.  If  9  is  the  angle  between  an  incident  ray  and  a  reflected  ray  then  cos(0)  >  —1. 

3.  There  are  no  rays  returning  to  the  surface  y  =  0  and  grazing  this  surface. 

Note  that  the  latter  two  conditions  in  the  SRG  hypothesis  ensure  according  to  Rakesh 
[12]  that  the  forward  map  T  is  a  Fourier  Integral  Operator. 

In  particular,  all  of  the  disks  fij  defined  in  section  3  are  assumed  to  lie  in  the  set 

n°(x„T). 

Our  next  goal  will  be  to  establish  the  following  result. 


Theorem  4.1  Assume  that  when  rays  are  traced  through  a  medium,  that  a  caustic  A 
forms.  Then  under  the  simple  ray  geometry  hypothesis,  the  forward  map,  T  is  bounded 
from  J3rc1/2(n°(x„T))  to  L2([0,L\  x  [0,T])  iff  0O°(xs,T)  fl  A  =  0. 


The  operator  T  is  in  fact  a  Fourier  Integral  Operator  (FIO)  and  hence  the  proof  of  the 
theorem  requires  a  knowledge  of  FIO’s.  An  introduction  to  this  theory  can  be  found  in 
the  following  references  [4], [17], [8], [6], 
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We  note  also  that  even  when  caustics  are  produced,  the  operator  T  (which  is  now 
unbounded)  still  retains  the  property  of  being  an  FIO  (  of  order  1/2  )  provided  that  there 
are  no  rays  leaving  the  source  and  returning  to  graze  the  plane  y  =  0.  This  latter  result 
is  shown  in  [12].  In  the  latter  case,  it  also  has  to  be  assumed  for  technical  reasons,  that 
the  source  location  x,  is  strictly  below  the  surface  i.e.  has  positive  y,  coordinate. 

A  large  part  of  this  section  will  be  devoted  to  the  proof  of  the  theorem.  We  will  need 
to  give  some  background  to  establish  notation.  We  will  also  need  several  lemmas  before 
we  commence  the  main  proof  of  this  theorem. 

First,  we  recall  an  important  fact  regarding  a  subclass  of  FIO’s  called  Pseudodifferential 
Operators  ($DO’s)  -  see  [4],  page  154.  Suppose  P  is  a  ¥DO  of  order  fi  then  it  has  the 
following  continuous  mapping  property  on  the  Sobolev  spaces 

P  :  HI  -  HZ?. 


It  will  be  useful  to  study  the  so  called  L 2  normal  operator  M  which  is  simply  defined 

by 

M=  *To  T 

where  XT  is  the  L2  adjoint  of  T . 

Denote  by  Fj,Afj  the  restriction  of  the  forward  and  normal  maps  respectively  to  the 
disks  fij  defined  in  section  3,  so  that  we  have 

J rj:^/a(fij)^i*([0,I]x[0,r]) 

'fj  :  L\[0,L}x[0,T})~ 

AT,  :  H]l\Uj)  ~  tfc-1/2(^).  (4-2) 

We  remark  that  we  will  only  use  6cj  £  Cq°  but  note  that  this  function  space  is  dense  in 
the  space  of  distributions  and  in  particular  dense  in  i/i^fij).  The  idea  being  that  if  the 
forward  map  were  to  be  discontinuous  (unbounded),  then  it  should  be  unbounded  on  the 
test  functions  too  (by  the  density).  In  fact  we  will  construct  a  sequence  of  test  function 
reflectivities,  and  then  show  that  the  forward  map  is  unbounded  on  this  sequence  in  the 
sense  stated  in  the  theorem. 

We  can  now  start  to  state  and  prove  the  lemmas  needed  to  prove  the  above  theorem. 
Lemma  (7.1)  below  gives  an  asymptotic  expansion  for  A/},  which  is  derived  in  [14]. 

For  a  source  location  x,  and  a  receiver  location  xr  on  y  =  0  define 
0(xs,  xr,  x)  =  r(x„  x)  +  r(x,  xr), 


12 


Trip  95 


V<f>  =  |V<^|  (  sin(V0)Cos(^0  ), 

V^€i?n,  i=  ||. 


(4.3) 

Also  given  the  pair  (x,£),  there  is  a  unique  location  on  y  =  0,  xr  =  xr(x,  f)  defined  as 
follows.  Consider  a  ray  leaving  the  source  x,,  arriving  at  x,  and  reflecting  so  that  —  ( is  the 
bisector  of  the  angle  0(x,,xr,x)  made  by  the  incident  ray  and  the  reflected  ray  emerging 
at  xr  -  see  figure  (4.2). 

Denote  by  6c j  the  Fourier  transform  of  6c j.  Then  we  have 


Lemma  4.1 


where 


A/jtfcj(x)  =  e'*(nj(xyZ)6cj(t), 

nj(x,0  =  n0,j(x,0  +  O(x,0 

(c~$(x)m(x„Xf,tf)a(x„x)a(Xr,x)) 

0,j(X7  0  l£l  '  (1+cosfl(x<1Xr,X))^^:(X.,Xr,X) 


Xp=Xr(X,0 


(4.4) 


Lemma  (4.1)  says  that  the  normal  operator  A/j  is  a  Pseudodifferential  operator  (4' DO) 
[15]  (of  order  1). 

The  following  result  will  become  important  to  us  soon. 


Lemma  4.2 
fij  such  that 


There  exists  a  conic  set  E  in  R2  —  {0}  based  at  0,  and  a  sequence  of  disks 

inf  -  oo  as  j  ->  oo  ,  V  {  €  E  n  S1.  (4.5) 

xen,  |{| 


Proof.  Define  i?j(x)  as  the  characteristic  function  of  the  disk  Slj.  Then  according  to 
Lemma  (4.1),  we  have 


rij, o(x,Q  _  (c  *(x)m(x<>xr,tr)a(x,,x)a(xr,x))  ^ 
If  I  (1 +  cos0(x„xr,x))?|£(xs,xr,x)  3 


(4.6) 


Xr=Xr(X,0 
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The  function  i?j  has  been  appended  to  emphasise  the  fact  that  it  is  understood,  x  €  ftj, 
where  flj  is  defined  as  in  section  3.  Towards  the  end  of  this  proof  we  will  modify  the 
definition  of  ft  j  slightly. 

The  amplitude  term  o(x,,x)  in  (4.6)  is  unbounded  for  an  incoming  wave  with  un¬ 
bounded  frequency  content  as  discussed  at  the  end  of  section  3.  The  same  is  true  for 
a(xr,x)  because  it  satisfies  the  transport  equation  (2.1)  which  can  be  written  as  an  O.D.E. 

x(l))+iv-Vr(x,,x(i))  =  0. 

Since  a(xr,x)  obeys  an  ordinary  differential  equation  (transport  equation)  with  a  nonzero 
initial  condition  we  can  bound  this  amplitude  from  below  by  some  constant. 


x2  x. 


/ 

Portion  of  a  Cauatic 


Figure  4.1:  Argument  for  boundedness  of 

With  regard  to  the  term  i?j(x)^-(x„xr,x),  we  must  be  sure  that  this  term  does  not 
become  large  with  increasing  j.  For  then  it  could  begin  to  compete  with  the  amplitude 
terms,  and  the  unboundedness  the  ratio  in  (4.6)  would  no  longer  be  obvious.  Therefore 
we  give  the  following  argument  to  prove  that  if  we  restrict  the  cone  in  which  £  varies 
then  we  can  always  avoid  this  situation.  Pick  a  point  xc  on  the  caustic  A  such  that  there 
is  a  ray  R\  starting  at  xc  and  intersecting  the  plane  y  =  0  transversally  at  x,  -  this  is 
possible  by  the  definition  of  a  caustic.  Denote  by  aj  the  takeoff  angle  from  the  vertical 
made  by  Ri.  Now  because  ft°(xs,T)  is  an  open  subset  of  y  >  0,  there  must  exist  a  point 
y  G  A  sufficiently  close  to  xc,  and  an  angle  02  sufficiently  close  to  01,  such  that  the  ray 
R2  departing  from  y  at  an  angle  of  a?  must  intersect  R\  at  some  point  z  €  ft°(x,,r) 
and  also  R2  stays  within  ft°(x,,T)  and  intersects  y  =  0  at  some  point  X2  transversally. 
The  latter  observation  is  true  because  of  continuous  (smooth)  dependence  of  solutions  to 
O.D.E.’s  (rays)  on  the  their  initial  data.  If  x2  =  xs  then  the  SRG  hypothesis  is  violated  at 
z.  Therefore  we  can  assume  that  the  ray  R2  emerges  on  y  =  0  at  X2  ^  x,.  Now  we  have 
the  following  situation  (see  figure  (4.1)  ):  there  are  two  rays  departing  from  y;  namely 
the  ray  R2  and  the  ray  R$  defined  by  the  ray  connecting  x,  to  y  (this  exists  since  y  is  a 
caustic  point).  Let  R3  have  takeoff  angle  0:3,  so  that  the  rays  R2  and  U3  depart  from  y 
with  take-off  angles  02,  <23  respectively,  and  arrive  back  on  y  =  0  at  two  different  locations. 
Therefore  there  must  exist  a  ray  departing  from  y  with  take-off  angle  a*  6  (012,013)  and 
intersecting  y  =  0  at  some  point  x*  such  that  ^^(x^xj^y)  ^  0.  In  fact  y  could  be  any 
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caustic  point,  so  ^(x,,x*,xc)  is  finite,  provided  the  (  varies  in  a  sufficiently  small  conic 
set  E  which  has  axis  equal  to  the  line  that  makes  an  angle  a*  with  the  vertical.  Therefore 
|£(x„xr(x,0,x)  is  bounded  independent  of  j  if  we  restrict  the  cone  in  which  £  varies, 
or  correspondingly  restrict  the  interval  in  which  xr  varies. 

Also,  for  any  j ,  we  can  always  pick  a  suitable  L  and  T  such  that  m(xr(x,  £),  tT)  =  1, 
whenever  x  €  ftj,  £  €  S  D  51  and  tr  €  [0,  T]. 

Clearly  all  other  terms  in  (4.6)  are  such  that  we  can  find  a  sequence  of  positive  numbers 
{Cj},withCj  — ►  oo,  so  that  for  (x,^)  €  fyxS,  then  no,j(x,()  /  |£|  >  Cj.  This  completes 
the  proof  of  the  lemma.  □ 

Now  define  a  new  operator  Afj,  which  is  “extension”  of  Afj  got  by  giving  it  the  same 
symbol  as  A fj  in  flj  x  E  and  extending  it  to  all  of  Qj  x  S1  so  that  its  principle  symbol 
no,j(x,£)  satisfies  the  inequality 


inf  nQtj(x,i)>Cj/ 2.  (4.7) 

xeo,, 

and  such  that  hj  is  still  a  symbol  of  order  one.  This  construction  is  trivial  to  carry  out,  and 
makes  Afj  a  strictly  elliptic  4* DO  of  order  zero.  Thus  if  Scj  has  Fourier  transform  rapidly 
decaying  outside  of  the  cone  E,  then  A fj6cj(x)  =  J\jScj(x)  modulo  a  smooth  function  of 
x. 


From  now  on  in  our  demonstration  of  unboundedness  of  the  forward  map  we  will  assume 
that  our  reflectivities  have  rapidly  decaying  Fourier  transform  outside  of  E  and  we  will 
accordingly  replace  A/j  by  Afj.  Under  this  assumption  Afj  differs  from  Afj  by  a  smoothing 
operator.  However  we  continue  to  write  Afj  in  place  of  Afj  for  convenience  of  notation. 

Proof  of  Theorem  (4-1): 

[=>]  We  first  prove  the  only  if  part  of  the  theorem.  To  prove  it  we  will  take  the 
contrapositive  point  of  view  and  assume  that  dil°(x,,T)f\X  ^  0.  Then  this  part  of  the 
theorem  will  be  proved  by  showing  that  T  is  unbounded  under  this  assumption. 

In  the  following  <  •,•  >  will  denote  the  usual  L2  inner  product  and  <  •,•  >(r,_r)  will 
denote  the  usual  Sobolev  space  duality  pairing.  Also  |j  •  ||a  will  denote  the  norm  in  the 
Sobolev  space  Ha . 

The  key  to  the  proof  of  the  theorem  is  Garding’s  Inequality  [15]  (page  55),  which  says 
that 

<  Afj6cj,Scj  ^  C\ J| Scj || d"  ^  Sj6cj,6cj 

where  Sj  is  a  smoothing  operator  that  arises  in  the  proof  of  the  Garding  Inequality.  This 
says  we  have 


II  Ti*cj 


<  6cj,AfjScj  >  >  CxJSc.Wl  -  Cr<j\\8cj\ 

2 


(4.8) 


15 


Trip  95 


We  will  apply  Garding’s  Inequality  to  Mj  (notice  that  our  error  of  replacing  Nj  by  Mj  could 
be  absorbed  by  the  smoothing  operator  Sj,  under  our  assumption  that  our  reflectivities 
have  Fourier  transform  decaying  rapidly  outside  of  the  cone  5). 


by 


The  proof  of  Garding’s  Inequality  says  that  r  can  be  any  real  number,  and  is  given 


'i»j 


1  inf 


(4.9) 


2  xenj.Ves1  l£l 

(we  have  just  shown  that  this  sequence  is  unbounded). 

Next  let  us  fix  r  =  0. 

Our  goal  now  is  to  show  that  we  can  construct  6c j  (by  choosing  \kj\  appropiately)  such 
that  the  first  member  of  the  right  hand  side  of  (4.8)  is  the  dominant  term. 


Let  e  >  0  and  for  any  j  €  N  we  want  to  find  6cj  such  that 

„  _  Iloilo, 


Cqj  <  t  ,  Vj. 


(4.10) 


Having  succeeded  in  showing  this  is  possible  then  the  first  term  on  the  right  hand  side  of 
(4.8)  is  dominant  and  we  get  unboundedness  of  the  quantity 


_  (6chMj6ci)  _  \\F6cjWl 

3~  M 


(4.11) 


I!  Ml 

2  2 

Now  we  need  just  prove  that  we  can  find  a  reflectivity  sequence  such  that  the  relation 
(4.10)  holds. 


Recall  that  our  reflectivity  sequence  has  the  form 

<MX)  =  ACi(x)e’fc,"x 

where  \j  €  Then  we  have  that 

IIM1  =  /  Ai(OI«i,(C)|!  K 


(4.12) 


where  A«)  =  (1  +  |C|2)* 


=  /  /  e'X  (kj  °Xj(x)dx 


dC 


=  /  A*(C  +  *i)|/  e-x  CXi(x)  dx| 
So,  =  J  A?(C  +  MIXi(C)|2 


dC. 

dC 
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Now  use  Peetre’s  Inequality  [13]  (page  21), 

As(u;)  <  2WAW(u>  -  ri) Xa{T}) 
with  s  =  |,  u  =  kj,  t]  =  kj  +  (  to  get 

>  2-Ui(tj)  J  A-t(C)  lxj«)l3<*C 


Therefore  we  have 


r;  <  C0j 


llXillli. 


—  C  01  \(le  llxjllo 

0J  (,)  ItallV 

2 

Next  we  choose  the  magnitude  of  kj  sufficiently  large  so  that 

W  1^-  Co,,  <  £,  Vj 

”2 

holds;  thus  yielding  rj  <  c  independently  of  j.  Note  this  relies  on  the  fact  that  the  ratio 
H^|jj0  is  independent  of  kj.  This  proves  the  “only  if”  part  of  the  theorem. 

—  J 


[<=]  The  “if”  part  of  the  theorem  is  trivial,  because  if  9ft0(x,,  T)  p|X  =  0,  we  know 
that  the  normal  operator  XTT  is  a  ¥D0  (see  [16],  [14]).  This  means  that  it  is  a  continuous 
operator  on  the  Sobolev  spaces,  as  remarked  before.  The  result  follows,  and  the  proof  of 
the  theorem  is  complete.  □ 

This  was  the  hardest  result  to  prove  in  this  section  and  now  we  can  improve  on  Theorem 
1  in  the  form  of  the  following  corollories. 


Corollary  4.1  Assume  that  when  rays  are  traced  through  a  medium  that  a  caustic  A 
forms.  Let  O  be  a  subset  of  the  halfspace  y  >  0.  Then  under  the  SRG  hypothesis  we  have 
that  the  forward  map  T\u  is  bounded  from  — *■  L2([ 0,  L)  X  [0,T])  iff  ilf )  A  =  0. 


N 

Proof.  In  general  ^|n  is  a  locally  finite  sum  of  FIO’s  .F|n  =  ^  J7,-,  where  each  J7,  has 

i=i 

canonical  relation  given  by  a  canonical  graph  (see  [16], [11]).  The  proof  now  follows  from 
Theorem  3.1  and  the  proof  of  Theorem  4.1.  □ 
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Corollary  4.2  Assume  that  when  rays  are  traced  through  a  medium  that  a  caustic  A 
forms.  Let  SI  be  an  open  subset  of  the  “post-caustic”  region  i.e.  the  multiply  covered  region 
of  the  medium.  Then  ^jn  is  bounded  from  Hl^2(Sl)  — *  £2([0,  Z]  x  [0,r])  iffSlf\X  =  0. 

N 

Proof.  As  in  the  previous  corollary  T\n  —  ^  Ti  holds.  The  proof  of  sufficiency  follows 

t=i 

from  the  previous  corllary.  For  the  necessity,  we  construct  a  reflectivity  sequence  6c j  — ► 
c  €  H^2(Sl)  as  in  theorem  4.1.  However  this  time  let  us  arrange  for  the  following  to  hold. 
At  any  point  x  €  SI  there  several  sheets  of  the  Lagrangian  A  manifold  associated  to  T 
lying  over  x  i.e.  7r-1(i)  =  {Ai,  •  •  -,AP}  where  each  A;  is  associated  to  a  different  sheet  of 
A.  Let  us  arrange  that  the  wavefront  set  (see  [9]  for  definition)  of  c  only  intersects  one 
sheet  of  the  Lagrangian  A: 

VFF(c)p|A  is  connected.  (4.13) 

This  sheet  corresponds  to  the  canonical  relation  of  the  FIO  T\  for  example.  The  proof 
now  follows  as  in  theorem  4.1  because 

Mjbcj  =  tTxT\6ci  +  SjScj  (4.14) 

where  Sj  is  a  smoothing  operator  which  can  be  absorbed  into  the  smoothing  operator 
appearing  in  the  application  of  Garding’s  inequality  in  theorem  4.1.  The  proof  is  completed 
by  applying  theorem  4.1.  to  the  operator  T\.  □ 

Remark:  Note  that  the  question  of  boundedness  as  in  the  above  theorem  and  corol¬ 
laries  when  a  reflectivity  field  has  its  support  only  on  a  caustic  is  still  an  open  (no  pun 
intended!)  question. 

Remark:  Note  that  in  theorem  4.1  and  the  above  corollaries  we  have  an  interesting 
result  when  caustics  form  in  a  medium.  One  might  be  concerned  about  how  likely  it  is 
that  a  caustic  will  form.  In  the  paper  by  White  [18],  it  is  argued  from  a  stochastic  point  of 
view,  that  even  in  a  mildly  heterogeneous  medium  a  caustic  will  form  (after  a  sufficiently 
long  time)  with  probability  one! 
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y  m  ®  Source 


Figure  4.2:  Snell’s  Law  and  the  definition  of  xr(x,£). 


5  Numerical  Tests 


There  is  also  another  forward  map  related  to  F .  It  is  the  map  resulting  from  the  more 
general  source  f(t)6(x  —  xs)  on  the  right  hand  side  of  (1.6).  We  call  this  forward  map  F\f\ 
to  emphasize  its  dependence  on  /.  Since  F  results  from  an  impulsive  source  S(t)6(x  —  x„), 
then  we  must  have  F[f]  —  F  *t  /(•)•  ^  *s  obvious  then  that  F[f]  is  unbounded  too 
for  a  suitable  source  function  f  £  L2.  The  choice  of  /  is  important  in  the  numerical 
experiments,  since  it  is  the  modified  map  F[f]  that  we  will  implement  in  this  section.  For 
convenience  of  notation,  we  will  write  F  in  place  of  F[f]  throughout  this  section. 

In  order  not  to  have  to  figure  out  what  reflectivity  fields  give  a  large  ratio  Rj  (defined 
by  (4.11)  with  F  replaced  by  F[f]),  we  decided  to  to  run  a  power  method,  so  that  the 
computer  works  this  out.  First  though,  we  need  to  clarify  how  the  power  method  applied 
in  our  case. 

Define 

Mf  =  F]  o  Fj, 

where  F*  refers  to  the  Hc{Qj)  Hilbert  space  adjoint  of  Fj. 

This  quantity  is  of  interest  because  of  the  following  relation 

ll-^^jllo  =<  ATjScjJcj  >i  . 


Therefore  we  have  that  the  largest  value  of  |{o  in  the  unit  sphere  of  H2($lj)  is 

equal  to  the  largest  Eigenvalue  of  Mj  (recall  that  computationally  Afj1  as  a  linear  operator 
is  some  finite  dimensional  matrix!) 
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It  is  this  Hilbert  space  normal  operator  that  we  apply  the  power  method  to.  So  next 
we  calculate  T’  in  terms  of  quantities  that  we  can  implement. 

Define  the  following  operator  on  Sobolev  spaces 

A  =  (/  —  V  •  V)a, 

where  I  is  the  appropiate  identity  operator,  for  the  Sobolev  space  in  question. 


Lemma  5.1 

Mf  =  A"1  o  xTi  o  Tj 


Proof:  Let  Scj  6  H?  (Qj),  Vj  €  L2([0,  L\  x  [0,  X1]).  Then  we  have 

(Jv,'6Cj,  Uj)0  =  ^Cj,  37}v}')^_^  (5-1) 

Note  that  the  right  hand  side  of  the  last  equation  is  just  a  bounded  linear  functional  on 

Hc(il)  and  so  by  the  Riesz  representation  theorem  we  are  guaranteed  the  existence  of  the 
operator  T*  such  that 

{T ,6cj, Vj)Q  =<  Scj^FjVj  >i  (5.2) 

Thus  we  get 

<  iehJ*Vj  >i=  (fcjfTjVj)  _ 

V  2  »  2  f 

=»  _k) 

=>  A r}  =  T]Tj  =  A-1  'TjT,.  (5.3) 

This  concludes  the  proof.  O 

Four  experiments  were  performed  and  in  each  case  the  following  was  the  procedure. 
The  (Hilbert  space)  normal  operator  was  composed  of  several  operators  reflecting  the 
result  of  the  last  lemma.  These  maps  in  turn  were  implemented  by  finite  difference  codes 
or  in  the  case  of  the  A  operator  this  was  done  by  using  discrete  Fourier  transforms.  This 
normal  map  was  run  many  times  with  the  output  of  the  normal  map  being  normalized  and 
then  used  as  the  input  again.  In  this  way  the  power  method  was  finding  the  Eigenvector 
(reflectivity)  corresponding  to  the  largest  Eigenvalue  of  Afj1  (corresponding  to  a  large  value 
of  Rj). 
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It  was  also  noticed  that  in  all  cases  that  as  we  moved  Slj  away  from  the  source  the  ratio 
Rj  decreased  until  it  started  to  approach  the  caustic  when  it  started  to  increase  in  the 
vicinity  thereof  again.  On  the  other  side  of  the  caustic  it  of  course  decreased  again.  It  was 
found  that  just  off  centre  from  the  peak  of  the  caustic  (see  figure  5.9)  that  Rj  was  larger 
than  at  the  peak  which  also  seems  interesting;  perhaps  some  interference  is  responsible  for 
this  effect.  This  is  the  location  which  we  chose  for  our  tests  (see  figures  (5.1), (5.2),  etc). 
Note  that  we  implemented  Afj,  which  is  Af  restricted  to  fly,  by  injecting  fly  into  a  larger 
grid,  applying  AT  and  then  projecting  the  result  onto  Qj. 

The  four  experiments  performed  were  the  following.  In  the  first  experiment  a  peak 
frequency  of  126  Hz  was  used  for  the  source  in  a  medium  which  has  a  lens  (low  velocity 
zone)  in  the  background  velocity  cq  which  forms  a  caustic  (see  figure  (5.9)).  The  reflectivity 
got  from  the  power  method  is  shown  in  figure  (5.1)  and  the  output  of  T  is  shown  in  figure 
(5.3).  The  second  experiment  is  identical  except  the  lens  is  removed  to  leave  a  constant 
background  of  1500m/s  (the  lens  had  a  minimum  velocity  of  1200m/s  approx.)  The  results 
are  shown  in  figures  (5.2)  and  (5.4).  The  third  experiment  is  the  same  as  the  first  except 
the  source  peak  frequency  is  reduced  to  63  Hz;  see  figures  (5.5)  and  (5.7)  for  the  results. 
The  fourth  experiment  is  the  corresponding  lower  frequency  experiment  as  the  second;  see 
figures  (5.6)  and  (5.8)  for  the  results. 

The  shape  of  the  reflectivity  got  in  these  experiments  makes  sense  because  it  is  oscil¬ 
lating  in  a  vertical  direction  which  is  the  direction  of  the  impinging  rays  for  the  most  part 
and  hence  it  will  reflect  most  of  the  energy  vertically  back  to  the  receivers.  The  output 
of  T  is  easily  interpreted  too;  we  see  that  it  has  several  “layers”  in  it  corresponding  to 
reflection  from  the  various  peaks  in  the  reflectivity. 

For  the  sources  we  used  Ricker  wavelets,  which  roughly  speaking  is  the  second  derivative 
of  a  Gaussian  function.  More  precisely,  if  upk  denotes  the  peak  frequency  of  the  wavelet, 

/(<)  =  ^  w  =  (™P*)2- 

See  figure  (5.10)  for  an  example  of  a  ricker  wavelet.  In  our  tests  the  wavelets  were  nor¬ 
malized  so  that  the  sup  norm  of  their  Fourier  Transform  was  unity. 

The  results  for  the  maximum  values  of  Rj  are  indicated  in  the  captions.  These  results 
indicate  a  stronger  reflection  in  the  vicinity  of  the  caustic  and  this  is  enhanced  at  a  higher 
frequency  as  expected  since  the  theory  was  based  upon  high  frequency  asymptotics.  The 
next  step  to  carry  out  is  to  refine  the  grid  so  that  we  can  have  smaller  reflectivities  which 
can  be  finely  sampled  to  oscillate  rapidly  as  seems  to  be  what  the  reflectivity  “wants 
to  do”  in  order  to  get  a  large  value  of  Rj.  The  unfortunate  aspect  of  this  is  that  the 
computational  time  goes  up  by  a  factor  of  around  8  every  time  we  halve  the  grid  size  and 
time  step.  Also  we  need  to  consider  different  sources  in  view  of  the  comment  at  the  end 
of  the  last  section. 

If  one  examines  the  asymptotic  behaviour  of  the  wave  field  theoretically  in  the  neigh¬ 
bourhood  which  actually  includes  a  piece  of  a  cusp  shaped  caustic,  as  was  done  in  [10] 
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one  finds  that  the  amplitude  is  larger  than  the  field  away  from  the  caustic  by  a  factor 
of  w1/6  if  our  probing  wave  is  monochromatic  of  frequency  u.  This  is  not  a  very  large 
exponent  and  thus  indicates  that  quite  a  large  frequency  may  be  needed  to  actually  see  the 
effect  of  our  theorem.  This  again  is  limited  by  computational  resources;  especially  memory 
considerations  will  be  important  if  we  are  to  sufficiently  sample  our  source  function. 

The  following  pictures  of  The  have  time  along  the  vertical  axis  (increasing  downward) 
while  the  xr  coordinate  is  the  horizontal  coordinate.  For  the  reflectivities,  the  vertical  axis 
is  the  “y”  spatial  axis,  and  the  horizontal  axis  is  the  “z”  spatial  axis. 
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Figure  5.1:  Reflectivity  6c  j  -  Pk  Freq.  126  Hz;  Lens  Medium;  Rj  —  0.10e-2.  Horizontal 
axis  =  x,  vertical  axis  =  y. 
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Figure  5.2:  Reflectivity  6c}  -  Pk  Freq.  126  Hz;  No  Lens;  Rj  =  0.21e-3.  Horizontal  axis  = 
x,  vertical  axis  =  y. 
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Figure  5.3:  Tjbcj  -  Pk  Freq.  126  Hz;  Lens.  Horizontal  axis  =  xr,  Vertical  axis  =  time  tT. 


Figure  5.4:  fjScj  -  Pk  Freq.  126  Hz;  No  Lens.  Horizontal  axis  =  xr,  Vertical  axis  =  time 
tv 
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Figure  5.7:  FjScj  -  Pk  Freq.  63  Hz;  Lens.  Horizontal  axis  =  xr,  Vertical  axis  =  time  tr. 
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Figure  5.8:  TjScj  -  Pk  Freq.  63  Hz;  No  Lens.  Horizontal  axis  =  xT ,  Vertical  axis  =  time 

tr. 
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Figure  5.10:  Ricker  wavelet  with  u  =  3.9e4.  Horizontal  axis  =  amplitude  in  meters, 
vertical  axis  =  time  in  milliseconds. 
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6  Conclusions 


We  have  shown  that  the  forward  map  associated  to  the  Linearized  Wave  equation  is 
discontinuous  in  the  prescence  of  caustics.  This  has  important  consequences  for  Linearized 
Inversion  techniques.  Any  local  methods,  for  example  gradient  based  algorithms,  will  have 
obvious  difficulties  because  of  the  loss  of  differentiability  and/or  continuity. 

As  a  consequence  of  the  unboundedness  of  T  we  have  that  the  normal  map  cannot 
always  be  a  Pseudodifferential  Operator  because  if  it  were  it  would  have  an  order  m 
associated  to  it  and  this  would  mean  that  it  would  map  —*  continuously 

but  as  we  have  just  seen  the  operator  grows  without  bound  as  we  let  Slj  become  closer  to 
the  caustic  (for  a  sufficiently  high  frequency  of  incident  waves). 

The  reader  will  notice  that  we  needed  to  use  a  reflectivity  field  with  indefinitely  high 
frequency  content  in  order  to  prove  unboundedness.  This  is  really  the  reason  for  the  title  of 
the  paper;  we  can  get  lower  estimates  on  the  forward  operator  bound  which  grow  without 
bound  as  we  let  Slj  approach  a  caustic  (for  a  sufficiently  high  frequency  content  in  the 
reflectivity  field).  This  means  that  we  will  see  an  “anomalous”  reflection,  meaning  that 
we  will  see  a  large  signal  relative  to  the  size  of  the  reflectivity  that  reflected  it  back  to  the 
surface. 

Note  that  the  results  in  this  paper,  were  for  a  two  dimensional  medium,  but  corre¬ 
sponding  results  hold  in  any  dimension. 

Finally  the  current  numerical  results  are  consistent  with  the  theory.  The  extent  to 
which  this  phenomenon  may  be  observed  is  still  undetermined  at  the  moment.  Further 
tests  will  be  carried  out  in  order  to  try  to  settle  this  question. 


7  Appendix 


Existence  of  Plane  wave  reflectors 

In  this  section  we  will  establish  the  intuitively  obvious  fact  that  the  plane  wave  reflectors 
mentioned  in  the  introduction  do  actually  exist  from  a  rigorous  point  of  view. 

Consider  the  setup  depicted  in  figure  (7.1).  First  we  set  down  some  notation.  Let  us 
denote  by  X  the  two  dimensional  vector  (z,  y )  as  indicated  in  the  diagram.  Let  V  be  a  C° 
vector  field  which  evaluates  the  unit  vector  parallel  to  the  velocity  vector  of  a  ray  leaving 
the  source  and  passing  through  the  point  of  evaluation  (see  diagram  for  example). 

Now  specify  the  direction  in  which  you  want  the  plane  wave  to  be  reflected  i.e.  specify 
a  vector  r  =  (rj,^);  this  choice  will  have  to  be  restricted  a  little  later  but  is  more  or  less 
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arbitrary  in  general.  Then  let  No  €  T(IOiVO)RJ  be  of  the  form 

No  =  N(X o)  =  N(x  o,  ito)  =  -7  1 . 1) 

(The  reason  of  this  choice  is  that  eventually  N  will  be  the  normal  vector  of  the  graph 
generated  by  y  =  V’C*)?  for  some  C1  function  i>.)  Also  choose  N  in  such  a  way  that  the 
following  equation  which  expresses  “Snell’s  Law  of  Reflection”  is  satisfied 

N0-(r  +  V(XQ))  =  0 


thus  we  have  that 

-  rj£ 0  _  Vi(.Yo)£o  +  ^2  +  V2(Xq)  =  0  (7-1) 

In  general  we  want  this  latter  equation  to  hold  with  Xo  replaced  by  X  and  £0  by  £(*)  i-e. 
suppose  we  can  have  some  reflecting  curve  given  by  the  graph  of  y  =  then  the  unit 
normal  vector  to  it  is  given  by 


1 

v/TTFW 


(-*(*),!) 


where  £(x)  =  Vi/>(i).  Therefore  our  task  is  complete  if  we  can  solve  (locally)  the  equation 

F(x,£)=  Vi(x,y(i))(  +  r2  +  Vi(x,y(x))  =  0  (7.2) 

(with  £  =  £(x))  since  then  all  we  have  to  do  is  find  the  function  if)  by  solving  the  ordinary 
differential  equation 

’  vy>(*)  =  -£(*) 

■ 

v>(zo)  =  yo 


Orientation 


Figure  7.1:  A  Plane  wave  reflector 


The  question  of  local  solvability  of  the  equation  (7.2)  is  answered  by  the  Implicit  Map¬ 
ping  Theorem  which  answers  “yes”  to  the  local  existence  of  a  C°  function  £(x)  satisfying 
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equation  (7.2)  provided  that  §jr(Ao>£o)  #  0  be-  rx  +  Vi(Xb)  ^  0.  Thus  the  choice  of  ri 
is  restricted  in  a  trivial  way.  This  completes  the  proof  of  local  existence  of  plane  wave 
reflectors. 
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